Comparison of different ladder models 
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Using density matrix renormalization group calculations, we compare results obtained for the t — J , 
one-band Hubbard and three-band Hubbard models of a two-leg CuO ladder. Spin and charge gaps, 
pair binding energies, and effective pair hoppings are calculated for a wide range of parameters. All 
three models have an insulating state at a filling corresponding to one hole per Cu site. For physically 
relevant parameters their spin gaps are similar in size but they exhibit quite different charge gaps. 
We find that the binding energy of a pair of doped holes is of the order of the undoped ladder spin 
gap for all three models. The main difference between the models is the size of the effective pair 
hopping, which is significantly larger in the three-band model with parameters appropriate for CuO 
materials than in the other two models. 
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PACS Numbers: 74.20.Mn, 71.10.Fd, 71.10.Pm 

The 2-leg CuO ladder materials have provided an in- 
teresting testing ground for ideas originally formulated to 
describe the 2D cuprates.El In particular, the undoned, 2- 
leg ladder material SrCu20^ exhibits a spinj^apH and 
the doped {SrCa)i4^Cu24:Oii ladder materialH can be- 
come superconducting under high pressure. The reduced 
dimensionality of this system has allowed for detailed nu- 
merical studies of both the t — J and the one-band Hub- 
bard models ota 2-leg ladder. In fact, it was a numerical 
Lanczos studyQ of a 2-leg t — J ladder which first sug- 
gested that the doped system might exhibit super|£on- 
ducting pairing. Since then, analytic calculationdju as 
well as density matrix renormalization group (DMRG) 
calculationsLl on long one-band Hubbard ladders have 
shown the (i2.2_j,2-like structure of the pairs and their 
power law correlations. Here we extend the DMRG ap- 
proach to study a three-band Hubbard model of a CuO 
2-leg ladder with the goal of comparing the t — J, one- 
band, and three-band Hubbard models. 

The three-band Hubbard model we have studied has 
a hole Hamiltonian which can formally be thought of in 
terms of the "3-leg" ladder shown in Fig. 1. Alternating 
sites on the top and bottom legs have Cu{dx^-y2) and 
0{px) orbitals, while the center leg has 0{py) orbitals 
on the sites bridging the Cu{d^2_y2) orbitals. With the 
orbital phase convention we have chosen, all the O — Cu 
hopping matrix elements arc — tp^. The difference in the 
hole site energies is A — Sp — Ed, and there is an onsite 
Cu Coulomb interaction Ud- The undoped system has 
one hole per Cu. In this state, if Ud is large compared to 
A, the charge gap is setJpy A and the system is said to be 
a charge gap ins.uIa.tor.u Various parameter values have 
been suggested,Lrtj with typical ones having A/tpd —2 
to 3, Ud/tpd — 8, and tpd ranging from l.SeV to 1.8 
eV. Throughout this paper we will confine ourselves to 
an isotropic ladder and leave the case in which the rung 
paranieters differ from the leg parameters to a future 
study.y 

The 2-leg Hubbard model has a hopping —t between 
the near neighbor ladder sites (along the legs and across 
the rungs) and an onsite Coulomb interaction U. Here 



the undoped state corresponds to half-filling. At half- 
fiUing, the charge gap, depends upon U arid is. said to be 
a Mott-Hubbard gap.l3 Typical parametersLltj are U/t = 
10 to 12 and t = 0.4 to 0.45eV. 

The t — J model has a near neighbor hopping term 
—t with a restriction that no site can have two fermions. 
The exchange interaction in the t — J model has the usual 
form 
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At half- 



filling, the t — J model is just the Heisenberg 2-leg laddp; 
and has an infinite charge gap. Typical parametersci t2l 
are J/t = 0.3 and t =0.4 to oJseV. 

We use DMRG techniquealJ to study long ladders (up 
to 64 X 2 sites in the t — J and one-band Hubbard mod- 
els and up to 16 X 2 Cu sites in the three-band Hub- 
bard model) with open boundary conditions. DMRG has 
been shostcutuiie a very accurate method to study ladder 



systems 
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lo Here we use up to 1200 states per block 




FIG. 1. Schematic of a CuO ladder. Here the solid circles 
represent Cu{d^2_y2) orbitals and the open circles represent 
0{px) orbitals along the upper and lower legs and 0{py) or- 
bitals on the rungs. There is a hopping matrix element —tpd 
between the O and Cu sites as shown by the solid lines. The 
energy difference between the O and Cu sites is A — Sp — Sd 
and there is an onsite Cu Coulomb interaction Ud- 
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and extrapolate DMRG results for energies to extract the 
limit of zero truiwation error and check the precision of 
our calculations.O In addition, a few tests performed on 
small clusters with periodic boundary conditions show, a 
perfect agreement with exact diagonalization results. ta 

We begin with the undoped ladder and calculate the 
spin gap 



A, = Eo {S, = l)-Ea (5, = 0) 



(2) 



Here Ag is the difference in energies between the spin 
1 andj-spin ground states. For the Heisenberg 2-leg 
laddeilia 
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As ~ - 
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(3) 



which is shown as the solid curve in Fig. 2, where we 
have plotted Ag versus At/ J . The dashed curve in Fig. 2 
shows results for a 32 x 2, half- filled Hubbard ladder. 
For large values oiU/t, the Hubbard model at half- filling 
maps to the Heisenberg model with J ~ At'^/U. Thus in 
Fig. 2, the two curves approach each other at large values 
of U/t. Using twice the spin gap as a measure of the 
effective exchange interaction, we see from Fig. 2 that 
the strong coupling expression J = At'^/U for the one- 
band Hubbard model overestimates the strength of the 
exchange interaction for physically relevant values of the 
parameters. However, we can use A^ as a unit of energy 
in comparing the t — J and one-band Hubbard models, 
each in their relevant physical parameter regimes. Thus, 
for an infinite 2-leg ladder, the t — J model with J/t = 0.3 
has Ag/i = 0.15 while the one-band Hubbard model with 
U/t = 12 has Ag/t ^ 0.12 as seen in Table 1. Taking 
t ~ 0.45eV, twice these spin gap energies give reason- 
able effective exchange couplings of 0.135eV (1600K) and 
O.lleV (1300K) for the t - J and one-band Hubbard lad- 
ders, respectively. Naturally, either of these two could 
be further adjusted by using a different value of t, but 
our point is simply that they are in the correct range. 
We also note that the spin gap A^/i = 0.11 of the one- 
band Hubbard model with U/t = 6 is similar to the value 
obtained for U/t = 12 (see Table 1). 
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FIG. 2. The spin gap As/t versus it/ J for the t—J (Heisen- 
berg) model (solid line) and U/t for the one-band Hubbard 
model (dashed line) at half-filling on a 32 x 2 ladder. 
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FIG. 3. (a) The spin gap As/tpd versus A/tpd for the 
3-band Hubbard model at a filling of one hole per Cu, with 
Ud/tpd = 8. The dashed line is the strong coupling limit, 
Eq. (^j. (b) The spin gap energy As versus Ud/tpd with 
A/tpd ~ 2 for the 3-band Hubbard model at a filling of one 
hole per Cu. All results are for an 8 x 2 Cu ladder. 

Results for the spin gap As of the three-band Hubbard 
model at a filling of one hole per Cu are shown in Figs 
3(a) and (b). In Fig. 3(a), Ag/tpd is plotted versus A/tpd 
for a ladder containing 8x2 Cu sites with Ud/tpd = 8. In 
strong coupling, the Cu — Cu exchange interaction has 
the form 



Jcu = 4 
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The dashed curve in Fig. 3(a) shows Jcu/ '2' and one sees 
that the spin gap of the three-band model approaches 
the strong coupling Jc«/2 result at large values of the 
A = Ep — £d splitting. However, just as for the one- 
band Hubbard model, the strong-coupling expression for 
the exchange interaction substantially overestimates it in 
the physically relevant region of parameters. The depen- 
dence of As on Ud/tpd for A/tpd = 2 is shown in Fig. 3(b). 
At large values of Ud/tpd, the exchange is dominated by 
the A term in Eq. (Q) as expected for a charge trans- 
fer insulator. Taking A/tpd = 2 and Ud/tpd = 8 we 
find that Ag/tpd extrapolates to = 0.035 for an infinite 
CuO ladder (see Table 1), which for tpd = 1.5eV gives 
Jcu = 2As = O.lleV (1300K). Thus, even though the 
region of relevant physical parameters does not lend it- 
self to a simple strong coupling expansion which relates 
the three models, the spin gaps for the three models are 
in fact similar in size and provide a useful energy scale 
describing the strength of the exchange coupling. 

The charge gaps for the insulating (undoped) phases, 
however, are quite different. The insulating state of the 
t — J model has an infinite charge gap reflecting the con- 
straint that no site can have two fermions. In the Hub- 
bard models, the charge gap is defined by 



2 



A, = {Eo{2) + Eo{-2) - 2Eo{Q)) /2, 



(5) 



where Eq (n) is the ground state energy of a ladder with n 
holes relative to the undoped ladder. The charge gaps for 
the half-filled, one-band Hubbard model and the three- 
band Hubbard model with one hole per Cu are plotted 
versus the onsite Coulomb interaction in Fig. 4(a). At 
large values of U, the charge gap in the one-band Hub- 
bard model varies as U, while as expected, the charge gap 
of the three-band Hubbard model saturates at a value 
set by A when Ud becomes very large. As noted above, 
the spin gap energy provides a useful unit of energy in 
comparing the different models. In Table 1 we show the 
charge gap in units of the spin gap in the one-band Hub- 
bard model for U/t = 6 and 12 and in the three-band 
Hubbard model for the physical parameters A/tpd — 2 
and 3 with Ud/tpd = 8. One can use either U/t = Q or 
12 to reproduce the three-band model Ac/Aj ratio de- 
pending on the value of A/tpd- Taking the same values 
of t and tpd as above, we obtain Ac = 1.4 to 3.9 eV. The. 
experimental value of the charge gap is still debatedtll 
but is of the same order of magnitude (2-5eV). 
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FIG. 4. (a) The charge gap Ac versus U/t for the one-band 
(dashed) Hubbard model on a 32 x 2 ladder and versus Ud/tpd 
for the three-band (solid) Hubbard model with A/tpd = 2 on 
a ladder containing 8x2 Cu sites. Note that Ac is measured 
relative to t for the one-band Hubbard model and relative to 
tpd for the three-band Hubbard model, (b) The hole occupa- 
tion per Cu (ncu) and O (no) versus the hole concentration 
1 + X per CuOi.5 unit cell on a ladder containing 16 x 2 Cu 
sites with Ud/tpd = 8 and A/tpd ~ 2. 



The "charge transfer gap" behavior of the three-band 
model is also seen in Fig. 4(b) where we have plotted 
the hole occupation per Cu and per 1.5 O versus the 
hole concentration 1 + x per CuOi.5 unit cell for the 
physical parameters Ud/tpd = 8 and A/tpd = 2. For 
a; < 0, the holes go primarily onto the Cu sites (72% 
of the hole density is on these sites) while for a: > 0, 
the O sites are favored with 88% of the additional hole 
density going to these sites. For example, for x = 0.25, 
the average hole density on the Cu site is increased by 
only 0.03 (4%) while that on an oxygen site is increased 
by about 0.15 (77%), or 0.22 per O1.5, from the result 
for the undoped ladder {x = 0). This behavior can easily 
be understood. For a; < 0, only Cu orbitals are occupied 
because they have lower energy than the O orbitals. (The 
hybridization of the Cu and O orbitals due to the finite 
hopping tpd is responsible for the fractional density on 
Cu and O sites.) For x > 0, all Cu orbitals are occupied 
by at least one hole. The energy to put a second hole 
on one of these orbitals is set by Ud, while the energy to 
put a hole on the O orbitals is set by A. Thus, for the 
parameters considered here (A < Ud), additional holes 
go onto O orbitals. In the opposite limit Ud < A, we 
have found that holes go primarily on Cu orbitals for 
x > 0, as expected. 

It is interesting to note that magnons also go primarily 
on Cu sites (76% of the spin density is on these sites). 
This distribution is very close to the hole distribution of 
the undoped ladder and remains constant for all doping 
studied —0.125 < x < 0.25. This result suggests that 
low-energy spin excitations involve only unpaired holes in 
(hybridized) Czi-arbitals even when x ^ 0. According to 
Zhang and Ricelia, each doped hole is locked in a singlet 
state with another holes in a Cu orbital and does not 
contribute to spin excitations at low energy. Our results 
certainly support this scenario although it is not obvious 
that their argument based on a CuO^ cluster is valid for 
the lattice configuration shown in Fig. 1. 

Next we consider the two-hole pair binding, defined by 



Apt ^ Eo{2) + Eo{0) ^ 2Eo{l) 



(6) 



if the quantity on the right-hand side is positive and 
Apb = otherwise. It should be noted that the depen- 
dence of the pair binding energy on system size is signifi- 
cant. Moreover, Apj, generally increases when the ladder 
length increases, while other quantities, such as the spin 
gap, decrease. Therefore, when comparing values of Apt, 
one should always keep the corresponding system size in 
mind. In Fig. 5(a) Apf, is plotted versus 4t/J for the 
t — J model and U/t for the one-band Hubbard model on 
a 32 X 2 ladder. Just as previously found for the spin gap, 
the pair binding energy for the t — J model approaches 
that of the one-band Hubbard model at large values of 
U/t. Furthermore, although the two models have very 
different charge gaps, the scale of their pair binding en- 
ergies in the physically relevant parameter region is set 
by As. This is clearly illustrated in Fig. 5(b), where we 
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FIG. 5. (a) The hole pair binding energy Apt/t versus At/ J 
for the t — J model (solid) and U /t for the one-band Hubbard 
model (dashed) on a 32 x 2 ladder, (b) The ratio Apt/ As in 
both models. 

show the pair binding energy in units of the spin gap as 
a function of U/t in the one-band Hubbard model and as 
a function of 4t/ J in the t — J model on a 32 x 2 ladder. 
We note that the ratio Apb/A^ is surprisingly similar 
in both models despite the widely different behavior of 
Aph and A_shown in Figs. 2 and 5(a). As previously 
discussed,Lrlij the pair wave function for both of these 
models has a dx2_y2-\ike form. Moreover, we have found 
that the structure of single holes and hole pairs in the 
one-band Hubbard model with UJ-L— 8 is identical to 
those observed in the t — J modelEj 

Results for the pair binding energy of the three-band 
Hubbard model are shown in Figs 6(a) and (b). For 
Ud/tpd = 8, the pair binding peaks for A/tpd = 2 and 
as shown in Fig. 6(b), the pair binding energy increases 
as Ud/tpd increases. Similar results for a^^Ilu^Os clus- 
ter were foundJrom Lanczos calculations. EJ'LiJ More re- 
cently, Martinc3 has discussed pairing on small clusters 
in terms of rehybridization associated with the charge- 
transfer channels. A recent Quantam Monte Carlo study 
of the three-band Hubbard modcEJ yields similar results 
for clusters with up to 6 x 6 Cu sites. For the physical 
range of parameters appropriate to the cuprates, one also 
has Aph ~ As in the three-band model. However, the ra- 
tio Apty/As extrapolated for a ladder of infinite length is 
clearly larger than in both other models with their typi- 
cal parameters, as shown in Table 1. A better agreement 
with the three-band Hubbard model results is obtained 



by taking U/t = 6 for the one-band Hubbard model. In 
addition, in the three-band Hubbard model the hole pair 
wave function has a (ij,2_j,2-like form as determined from 
measurements of the rung-rung and rung-leg pair field 
correlations. 

We have also measured the pair binding energy versus 
the hole concentration x relative to the undoped ladder. 
In this case the pair binding energy is defined by 

Apb{x) = Eo{2n) + Eo{2n - 2) - 2En{2n - 1), ti > (6') 

Apb{x) = Eoi2n) + Eoi2n + 2) - 2Eo{2n + 1), 7i< (6") 

Here x is equal to 2n divided by the number of holes in 
the undoped ladder. Results for the pair binding energy 
Apb relative to the undoped spin gap A^ are shown in 
Fig. 7(a) for the three models. There is no striking dif- 
ference between the three models for a; > 0. We note, 
however, that while the Hubbard model is particle-hole 
symmetric about half-filling, one can only add holes to 
the t — J model and we find no evidence for pair binding 
for X < in the 3-band Hubbard model with Ud/tpd = 8 
and A/tpd = 2. We are not sure of the exact nature 
of the ground state of the three-band Hubbard model 
with less than one hole per Cu site. In fact, we have 
not been able to investigate this regime as thoroughly as 
the X > regime because density matrix renormalization 
group calculations are much harder and less accurate in 
this case. However, our numerical simulations strongly 
suggest that the a; < and x > regimes are quite dif- 
ferent. As noted previously, holes go primarily on Cu 
orbitals for x < while doped holes are on O orbitals for 
X > (see Fig. 4(b)). Thus, it is possible that the ef- 
fective interactions between doped particles are different 
for both regimes. 
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FIG. 6. (a) The pair binding energy Apt versus A/tpd with 
Ud/tpd = 8 for the 3-band Hubbard model, (b) The pair 
binding energy versus Ud/tpd with A/tpd ~ 2 for the 3-band 
Hubbard model. These calculations are for an 8 x 2 Cu ladder. 
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TABLE I. Spin gap (in units of the bare hopping term t or tpd), and charge gap A^, pair binding energy Apt, effective 
pair hopping t^ff and effective magnon hopping Vef / (in units of the spin gap) obtained by extrapolating to a ladder of infinite 
length. 



Model 


Parameters 


A, 


Ac/A, 


Apfc/As 


W//As 


Vef f/ As 


t-J 


J/t = 0.3 


0.151 


oo 


0.71 


2.4 


5.0 


1-band 


U/t = 12 


0.116 


70 


0.83 


5.4 


8.5 


1-band 


U/t = 6 


0.111 


25 


1.4 


12 


24 


3-band 


A/tpd = 3, Ud/tpd = 8 


0.030 


58 


1.2 


12 


7.5 


3-band 


A/tpd = 2, Ud/tpd = 8 


0.035 


29 


1.6 


18 


11 



The decrease in the pair binding energy with x cor- 
relates with the decrease in the near-neighbor spin-spin 
correlations {Si ■ Sj) shown in Fig. 7(b). (Note that for 
the three-band Hubbard model we show the spin-spin 
correlations between near-neighbor Cu sites.) Thus, the 
decrease in the pair binding energy with increased hole 
concentration reflects the destruction of the underlying 
exchange correlations by the added holes. It should be 
noted that the differences between the values of (Si-Sj) in 
the three models are mostly due to the increase of charge 
fluctuations in going from the t — J model to the one- 
band Hubbard model and then to the three-band Hub- 
bard model. These charge fluctuations reduce the local 
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FIG. 7. (a) The ratio of the pair binding energy to the 
undoped spin gap versus hole doping. The diamonds are for 
a 32 X 2 f — J ladder with J/t = 0.3. The circles are for a 
one-band 32 x 2 Hubbard ladder with U/t — 12. The squares 
are for the three-band Hubbard model with Ud/tpd ~ 8 and 
A/tpd = 2 on a 16 X 2 Cu ladder, (b) The near-neighbor Cu 
spin-spin correlation function {Si ■ Sj) versus x for the three 
models. 



magnetic moments (5'^^) and thus the absolute value of 
(Si ■ Sj). However, as noted previously, the effective 
exchange coupling between spins is similar in the three 
models. 

We have also calculated the effective hopping param- 
eter teff of a hole pair from the dependence of the the 
pair energy obtained in ladders of different lengths. We 
define the pair energy by 



£;o(2)-So(0). 



(7) 



In ladders with open boundary conditions, the energy of 
the pair varies as 



L. 



(8) 



where the effective system length Lgj f differs from the 
actual ladder length L because of end effects. For a given 
system, L^f f can be determined from the wavelength A = 
2{Lef / + 1) of the charge density distribution of the pair. 
Fig. 8 shows a plot of the ground state energy of a hole 
pair versus (ie// + 1)^^ for the three different models. 
teff is equal to the slope of the hues in Fig. 8 divided 
by TT^. We have found that the difference 5L = Leff — L 
tends to a constant for large systems. Thus, in practice 
we can substitute L + SL for L^f / in Eq. ^ and use SL as 
a fit parameter. 
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FIG. 8. Plot of the hole pair energy versus {Leff + 1)~^ 
for the three different models. The notation is the same as in 
Fig. 1^ Actual ladder lengths L range from 5 to 16 Cu sites 
in the three-band Hubbard model and from 16 to 48 for the 
two other models. For each model, the zero of the energy has 
been set to the extrapolated value of the pair energy. 
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In Table 1, we list teff normalized with respect to As 
for the t — J, and the one- and three-band Hubbard mod- 
els with J/t = 0.3, U/t = 6 and 12, and Ud/tpd = 8, 
^/tpd — 2 and 3, respectively. We have also listed the 
pair binding energy Apt, in units of the spin gap energy 
Ag. Now, as previously discussed, the spin gap in the 
insulating case is of order J/2 and it sets the scale of 
the pair binding energy so that Ag and Api, for the three 
models are quite similar. However, the pair dispersion is 
enhanced in going from the t — J model to the one-band 
Hubbard model and further enhanced for the three-band 
Hubbard model if one uses parameters appropriate for 
CuO materials. We believe that this enhancement is as- 
sociated with the additional charge fluctuations which 
the one and three-band Hubbard models allow. The large 
enhancement in the three-band model reflects the fact 
that its charge gap is set by A when Ud is large rather 
than Ud- In Fig. 9 we show the effective pair hopping as 
a function of U/t in the one-band Hubbard model and 
as a function of 4t/ J in the t — J model. In the t — J 
model it seems that the effective pair hopping is always 
small {teff/t < 0.4). In the one-band Hubbard model, 
however, the pair dispersion is strongly enhanced as U/t 
decreases and for U/t = 6 the ratio t^ff/As is similar 
to the value obtained in the three-band Hubbard model 
with A/tpd = 3 and Ud/tpd = 8 (see Table 1). For large 
values of U/t, the pair dispersion of the one-band Hub- 
bard model approaches that of the t — J model as shown 
in Fig. 9. However, the difference between the pair dis- 
persion in both models seems to be of the order 4t/U 
when [/ — > oo. This difference could be due to next- 
nearest- neighbor hopping terms of the order At/U which 
are neglected in the derivation of the t — J model from 
the strong-coupling limit of the Hubbard model. 

An effective magnon hopping Veff can be calculated 
from the dependence of the spin gap on the ladder length. 
We use Eq. ^jwith Veff and As substituted for ig// and 
tp, respectively. In this case the effective ladder length 
is determined from the spin density profile. In Table 1 
Vef f normalized with respect to As is listed for undoped 
ladders. 
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FIG. 9. The effective hopping parameter for two holes is 
shown by the circles for the one-band Hubbard model and by 
the squares for the t — J model. 



The values of Veff/As obtained in the three models are 
of the same order of magnitude for the physically rel- 
evant parameters. We note that in this case the value 
obtained in the one-band Hubbard model with U/t — 6 
is much larger than the other results. In undoped ladders 
the effective magnon hopping obtained for the one-band 
Hubbard model approaches the result found in the t ~ J 
(Heisenberg) model, Veff = 2.5 J, at large values of U/t 
as expected. A recentj^erturbation calculation of two- 
leg Heisenberg laddersEj gives Ug// = 2, in satisfactory 
agreement with our numerical result. 

To summarize, all three models of course have a spin 
gap in the undoped phase. Furthermore, the parame- 
ters for the different models can be chosen to make these 
spin gaps comparable and in the correct physical regime. 
As noted, however, these three models have very differ- 
ent charge gaps with the t — J model having an infinite 
charge gap, the one-band Hubbard model having a Mott- 
Hubbard gap set by U and the three-band CuO model 
having a charge transfer gap set by A = £p — e^; for the 
physical parameter range of interest. The binding en- 
ergy of two added holes is basically set by the spin gap. 
However, Apt/ As = 0.7 to 0.8 for the t — J and one- 
band Hubbard ladder while Apb/Ag = 1.2 to 1.6 for the 
three-band Hubbard ladder with parameters appropri- 
ate for CuO materials. Thus the hole pairs are in fact 
bound more tightly in units of the spin gap energy in 
the three-band Hubbard ladder. The main difference be- 
tween the models is the size of the effective pair hopping 
which is significantly larger in the three-band model for 
the physical parameters. It is often assumed that these 
three models with parameters and band-fillings appro- 
priate for CuO materials describe the same low-energy 
physics. Our study has not revealed any fact which ex- 
plicitly contradicts this point of view. However, our re- 
sults show that neither the one-band Hubbard model nor 
the t — J model can reproduce all the three-band Hub- 
bard model results with a single set of effective param- 
eters. For instance, it seems that the pairing properties 
are better reproduced by the one-band Hubbard model 
with U/t = 6 than with the usual parameters U/t — 12. 

Finally, we believe that the larger dispersion signified 
by the effective pair hopping, which we found for the 
three-band Hubbard model, points to an important fea- 
ture of the physics contained in the charge transfer insu- 
lator. Pairs are less likely to localize in the three-band 
Hubbard model than in the other models for the physical 
parameters. This could have a significant effect in a two- 
dimensional lattice. For instance, in the two-dimensional 
t — J model doped holes tend to form ordered arrays 
called stripes, which seem to suppress superconductiv- 
ity, at least when the stripes are static£j A study of 
a three-leg ladder using both density matrix renormal- 
ization group methods and quantum Monte Carlo simu- 
lations has shown that similar stripes also appear in the 
one-band-Hubbard model for U/t > 6, but not for weaker 
coupling Jl3 Thus, it seems that the sharp increase of pair 
mobility observed for decreasing U/t in the two- leg lad- 
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der correlates with a transition from a striped ground 
state to a ground state without stripes in wider ladders. 
Therefore, we think that the larger charge fluctuations 
which the three-band Hubbard model allows (compared 
to the t — J model) should lead in the two-dimensional 
lattice to a reduced tendency for domain walls to lock up 
in static arrays with suppressed superconductivity. 
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